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Abstract 

Formic acid is known to act as a reduction agent for copper oxide. Its 
thermal uni-molecular decomposition was studied by means of DFT with 
special attention to reaction paths and kinetics. 
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1 Introduction 



Thermal decomposition of formic acid proceeds into the two possible products, 
CO + H2O or CO2 + H2, of which CO is known to be responsible for reducing 
behaviour. 

Columbia and Thiel provide an overview about the interaction of formic acid 
with transition metal surfaces pQ, including the respective reaction paths but 
nothing about reaction kinetics. The thermal uni-molecular decomposition re- 
actions have already been studied by Ko Saito et al. [2[. Their work includes 
ab-initio calculations at the Hartree-Fock (HF) and M0ller-Plesset Perturbation 
Theory (MP-2) level and a quite exhaustive list of referenced experimental work 
as well as their own experiments of Ar diluted HCOOH in a shock tube system. 
DFT calculations were previously performed by Jan W. Andzelm et al. [3J, 
CCSDT-1 calculations by Goddard et al. 0]. 

This work aims at simulating the reaction pathway of the decomposition of 
formic acid in the gas phase in order to learn about reaction energies and kinet- 
ics, but also to make predictions about reasonable process parameters for CuO 
reduction, and possibly to support the results of earlier work. It is meant as a 
starting point for a more thorough investigation of the reaction mechanisms of 
the reduction itself. 

2 Simulation Details 

The calculations of reaction energies were run with DMol 3 [5] in Materials Studio 
and Turbomole [6] . Reaction paths were approximated with Synchronous Tran- 
sit methods [7] as implemented in DMol 3 , the guessed transition states refined 
and verified using Eigenvector Following methods as implemented in DMol 3 and 
Turbomole. The parameters for the DFT calculations were set as follows. For 
DMol 3 the DNP basis set, PBE 8 functional and the default 'fine' settings 
were employed, i.e. SCF cycles were converged to an (RMS) energy change of 
1 • 10 _6 eV, geometries optimized to a maximum force of 2 • 10 _3 Ha/A and a 
maximum displacement of 5 • 10 -3 A. In Turbomole, the TZP basis set, b3-lyp 
[TOj HU H21 [T3j [H] hybrid functional and equal convergence parameters were 
used. 

3 Reaction Energies 

SCF Energies of the examined structures with corresponding zero point vibra- 
tional energies (ZPVE) are displayed in Figure [U activation energies are listed 
in Table [1] See Figure Q] for an illustration of the various energy and enthalpy 
terms. 

The reaction barriers for the two channels are nearly equal (DMol 3 predicts a 
slightly lower barrier for path 2, Turbomole a slightly higher one, but the dif- 
ferences are in the order of expected errors). 

Previous ab-initio calculations and measurements are discordant as to whether 
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Table 1: Energy differences between cis- and trans-HCOOH and their corre- 
sponding TS in kcal mol -1 




Reactand Transition State 



Figure 1: Definition of the various Energy and Enthalpy Terms 



the reaction barriers are similar (Hsu et al. (1982) [15], Andzelm et al. (1995) 
[3] and Goddard et al. (4j) or significantly different (Saito et al. (1984) [2] and 
Blake et al. (1971) [16]). The latter two predict a significantly higher (about 20 
kcal mol~l) barrier for path 2, and make this responsible for the higher conver- 
sion to CO + H2O. However, looking at the current data it seems the difference 
between the reaction paths is rather due to thermodynamic factors, see section 
[4] for explanation. 
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Figure 2: DFT energies of the investigated structures. Literature values (black) 
from Ko Saito et al. [2] do not include ZPVE. The dashed line illustrates a 
possible reaction path, the parabolas illustrate a harmonic approximation of 
vibrational energies, their slope is arbitrary. 
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4 Kinetics 



4.1 Theory 

Reaction kinetics are evaluated in terms of transition state theory, which 
is expected to give a good approximation in the case of thermal equilibrium 
(i.e. fast energy exchange of a molecule with the surrounding system). A more 
complex theory like RRKM is not employed because the large uncertainties in 
DFT calculations would still render the results inaccurate. It must be noted 
here, that transition state theory gives an upper bound to the rate constants. 
The equilibrium constant for the first part of the reaction K$ (from formic acid 
to the transition state) is 

t = C(jf*) = Z X —AEq/RT 

c(HCOOH) zhcooh 

where , chcooh are the concentrations of the transition state (TS) and formic 
acid, z x , zhcooh their partition sums, respectively, and AE is the difference 
of ground state energies between reactant and transition state. R is the ideal 
gas constant. z x may now be written as 

z + = 



I _ e -hv/k B r J 



where v is the frequency of the vibrational mode of the TS along the reaction 
coordinate and z\ is the partition sum without that vibration. Since the corre- 
sponding force constant is very low, the frequency will also be very low and the 
exponential function is developed into a series and all but the linear term are 
neglected 

t- 1 1 k B T 

Z Z h- e-Wfe B T ~ z h-(l- hv/k B T) Z% hv 

The TS has 2v chances per unit time to dissociate (y in each direction) . Due to 
the very low force constant, one may assume, that every chance for dissociation 
will be used. However, only half of the chances will lead to the reaction products, 
the others will lead back to the reactants. Thus, the rate constant is 

k = 2v-yK t 

= V . Z% .*BL 1 e -AE /RT 

hv zhcooh 



k B T z t 



. e -AE /RT 



h zhcooh 

= ksT . e A S(T)/R _ e -AH vib (T)/RT . ^AEscf/RT 

h 

which is known as the Eyring Equation. Here, S is the entropy, H vi b the enthalpy 
of the system due to vibrations (including ZPVE) and Esc f the approximated 
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ground state energy of the system. The temperature dependence of AS and 
AH is displayed in figures [3] and |4j AH v n, is defined as 

AH mb = AH- AEscf 
4.2 Results and Discussion 

Reaction enthalpies, vibrational spectra and molecular structures were found 
to agree well between the two programs, therefore all of the presented kinetic 
results are from DMol 3 output. 

The reaction rate constants are equal at about 280 °C , at lower temperatures 
the decarboxylation reaction is preferred, at higher temperatures the dehydra- 
tion reaction is. 

Thermodynamic contibutions have a dramatic effect on the reaction rates. Ne- 
glecting all thermodynamic quantities but ZPVE in the Eyring equation yields 
a much higher (1-2 orders of magnitude) rate constant for path 2 compared to 
path 1. The following figures (J3J 2] and [5]) display the temperature dependence 
of the two reaction enthalpy and entropy barriers and the rate coefficients. 
Reaction barriers are consistent with other ab-initio works, the reaction rates 
are consistent with experimental results from Blake et al. but several orders of 
magnitude smaller than the results from shock-wave experiments. 
In order to exemplify these results we apply them to reactor conditions as de- 
scribed by Thomas Wachtler [T7]. The reduction is carried out at 388K (115 °C) 
and 1.3mbar with a flow rate of 70 mg per minute for formic acid and 180 mg 
per minute for Ar as a carrier gas the formic acid gas will have a molar con- 
centration of about 1 ■ 10~ 2 ^. The rate constant for the dehydration reaction 
at 115 °C is 2.42 • 10 _21 s _1 . Under these conditions the formic acid has about 
0.4 seconds to decompose before it reaches the sample (assuming one litre of gas 
volume between the supply and the sample). After that time the CO concen- 
tration in the gas will be 9.86 • 10~ 24 A copper oxide film that consists 
half of CU2O and half of CuO and measures 10 nm • 1 cm ■ 1 cm contains about 
6 • 10 -8 mol oxygen. Therefore it would take 4 ■ 10 16 min just to flow by the 
stoichiometrically equivalent amount of CO. See Appendix [A] for the details of 
this calculation. 

Clearly, the thermal decomposition cannot be the only reaction mechanism for 
the reduction. 



6 



-3,6 




-4,8 -| , , , , , , , , , ,-L 

200 400 600 800 1000 

Temperature [K] mns-HcooHn.-rsi 

DMol 3 



Figure 3: Enthalpy and Entropy differences between trans-HCOOH and TS1 
function of temperature 
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Figure 4: Enthalpy and Entropy differences between cis-HCOOH and TS2 
function of temperature 
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Figure 5: Rate coefficients as function of temperature. Dashed lines are values 
from Blake et al. 1971 



Development of the concentrations over time at 750K 
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Figure 6: Concentrations of HCOOH, CO and CO2 over time at 750K 
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A Applying results to reactor conditions 

In this section we want to apply the results to actual reactor conditions and 
calculate how long it might take to reduce a copper oxide film. 
A copper oxide film that measures 1 nm • 1 cm • 1 cm has a volume of 10~ 6 cm 3 . 
If half of its volume is CuO and half of it is Cu 2 then the respective molar 
amount of oxygen (atoms) is 

CuO 5 • l(T 7 cm 3 • • — ^ = 4 • l(T8mol 

cm 3 80^ 

Cu 2 5 • 10~ 7 cm 3 • • - A } = = 2 • l(T8mol 

cm 3 143_E_ 

total 6 • l(T 8 mol 

The process runs at 1.3mbar and 115 °C (388 K). Formic acid is supplied with 
a flow rate of 70 with an Argon carrier gas flow of 100 seem (standard cubic 
cm, approximately 180 ^^)- 

We use the ideal gas equation to get the total gas concentration: 

n v niol 
P V = nRT c total = - = JL = 0.0403 — 

V Hi m? 

We compute the partial molar fluxes and the concentration of HCOOH in the 

reactor, 

70 mo l 

HCOOH = 1.5207 • 10" 3 

46.03^ min 

180^ , mol 

mol 

total 6.027 -lO" 3 — 
mm 

chcooh = . 07 • 1-5207 = 0.0101688^ 

4.506+ 1.5207 m j 

and the gas flow speed Q. 

6 Q27 mmo1 m 3 / 

Q tot = SHL-HCOOH = 0.1495 = 2.492- 

0.0403^ min s 



Thus, it takes about 0.4 seconds to exchange a litre of the gas. Assuming the 
gas volume between supply and the sample is just one litre, we can calculate 
the CO concentration in the gas when it reaches the sample: 

ceo = chcooh ■ (1 - e" fcAT ) 

= 0.0101688^(1 - exp{-2A ■ 10~ 21 i • 0.4s)) 

m A s 

= 9 .86.10- 2 ^ 
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the CO molar flux is then 

Qco = ceo ■ Qtot = 9.86 • lCr 24 ^ ■ 0.15— = 1.4786— 

m J mm mm 

Thus, the time until 6 • 10~ 8 moles of CO have reached the sample will be about 

6 • 10~ 8 mol 



1.48 • 10- 24 ^2l 



4 • 10 16 min 
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Development of the concentrations over time 
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